Transcriptome Sequencing on Treatment of Kidney Deficiency and Blood Stasis-thin Endometrium in Rats with Bushen Huoxue Prescription Through Enema / 中国实验方剂学杂志

ObjectiveTo explore the mechanism of Bushen Huoxue enema in treating the rat model of kidney deficiency and blood stasis-thin endometrium （KDBS-TE） by transcriptome sequencing. MethodThe rat model of KDBS-TE was established by administration of tripterygium polyglycosides tablets combined with subcutaneous injection of adrenaline. The pathological changes of rat endometrium in each group were then observed. Three uterine tissue specimens from each of the blank group， model group， and Bushen Huoxue enema group were randomly selected for transcriptome sequencing. The differentially expressed circRNAs， lncRNAs， and miRNAs were screened， and the disease-related specific competitive endogenous RNA （ceRNA） regulatory network was constructed. Furthermore， the gene ontology （GO） functional annotation and the Kyoto Encyclopedia of Genes and Genomes （KEGG） pathway enrichment were performed for the mRNAs in the network. ResultCompared with the blank group， the model group showed endometrial dysplasia， decreased endometrial thickness and endometrial/total uterine wall thickness ratio （P<0.01）， and differential expression of 18 circRNAs， 410 lncRNAs， and 7 miRNAs. Compared with the model group， the enema and estradiol valerate groups showed improved endometrial morphology and increased endometrial thickness and ratio of endometrial to total uterine wall thickness （P<0.05）. In addition， 21 circRNAs， 518 lncRNAs， and 17 miRNAs were differentially expressed in the enema group. The disease-related specific circRNA-miRNA-mRNA regulatory network composed of 629 nodes and 664 edges contained 2 circRNAs， 34 miRNAs， and 593 mRNAs. The lncRNA-miRNA-mRNA regulatory network composed of 180 nodes and 212 edges contained 5 lncRNAs， 10 miRNAs， and 164 mRNAs. The mNRAs were mainly enriched in Hippo signaling pathway， autophagy-animal， axon guidance， etc. ConclusionBushen Huoxue enema can treat KDBS-TE in rats by regulating specific circRNAs， lncRNAs， and miRNAs in the uterus and the ceRNA network.

Exploring the potential mechanism of artesunate in intervening with U87 cells and U251 cells with distinct therapeutic effects on the basis of transcriptome sequencing and network pharmacology / 药学学报

Artesunate possesses the potential of intervening with glioma, however, its pharmacological mechanisms remain unclarified. Firstly, the effects of artesunate on cell activity, proliferation and apoptosis of U87 and U251 human glioma cells were explored. It was found that artesunate exerted stronger inhibitory effects on the activity and proliferation of U87 cells than U251 cells. It could significantly promote apoptosis in U87 cells (P < 0.05), while only high dose of artesunate can promote that of U251 cells (P < 0.01), detected by Hoechst and TUNEL cell apoptosis staining. Further, the differential expression gene sets between artesunate-sensitive and non-sensitive cell line, as well the therapeutic effects-related genes of artesunate were obtained through transcriptome sequencing and differential data analysis by using the lysates of U87 and U251 cells before and after artesunate treatment, aiming to explore the molecular mechanism of distinct artesunate sensitivity to two types of cells. Then, key putative targets that related to therapeutic effects were screened by constructing the interaction network of differential genes of three above comparison groups, and calculating their topological characteristics. Pathway enrichment analysis showed that those key putative targets were significantly enriched in several signaling pathways that were closely associated with the main pathological changes of glioma, among which apoptosis-related activating transcription factor 4 (ATF4)-DNA damage induced transcript 3 (DDIT3)- polyadenosine diphosphate ribose polymerase 1 (PARP1) signaling axis was the most enriched in. Molecular docking indicated that artesunate had fine binding affinities with ATF4 and DDIT3. Above all, this study preliminarily revealed that ATF4-DDIT3-PARP1 signaling axis is the target pathway of artesunate intervening with U87 glioma cells, and PARP1 may be an important gene for U251 cells to develop resistance to artesunate. Our results not only provide fundamental experimental evidence for artesunate as a potential therapeutic drug in glioma treatment, but shed light into overcoming drug resistance in its clinical therapy.

Maackiain inhibits proliferation and promotes apoptosis of nasopharyngeal carcinoma cells by inhibiting the MAPK/Ras signaling pathway / 中国天然药物

Nasopharyngeal carcinoma (NPC) is the third most common malignancy with a high recurrence and metastasis rate in South China. Natural compounds extracted from traditional Chinese herbal medicines have been developed and utilized for the treatment of a variety of cancers with modest properties and slight side effects. Maackiain (MA) is a type of flavonoid that was first isolated from leguminous plants, and it has been reported to relieve various nervous system disorders and exert anti-allergic as well as anti-inflammatory effects. In this study, we demonstrated that MA inhibited proliferation, arrested cell cycle and induced apoptosis in nasopharyngeal carcinoma CNE1 and CNE2 cells in vitro and in vivo. The expression of the related proteins associated with these processes were consistent with the above effects. Moreover, transcriptome sequencing and subsequent Western blot experiments revealed that inhibition of the MAPK/Ras pathway may be responsible to the anti-tumor effect of MA on NPC cells. Therefore, the effects of MA and an activator of this pathway, tertiary butylhydroquinone (TBHQ), alone or combination, were investigated. The results showed TBHQ neutralized the inhibitory effects of MA. These data suggest that MA exerts its anti-tumor effect by inhibiting the MAPK/Ras signaling pathway and it has the potential to become a treatment for patients with NPC.

Analysis of intercellular communication patterns in facioscapulohumeral muscular dystrophy based on single-cell nuclear transcriptome sequencing / 国际生物医学工程杂志

Objective:To analyze the patterns of intercellular communication in facioscapulohumeral muscular dystrophy (FSHD) by single-cell nuclear transcriptome sequencing.Methods:Bilateral asymmetrical lesions mouth orbicular muscle of two patients with FSHD and mouth orbicular muscle of two healthy patients were selected. Six samples were obtained, and were divided into control group, mild group and severe group. The normal orbicularis muscle sample was collected from 2 healthy individuals (the control group). The muscle samples in the mild group were from two patients with relatively normal muscle sides, and the samples in the severe group were from two patients with more severe muscle damage sides. Single-cell nuclear transcriptome sequencing was performed on all cells of the three groups. Reduced dimension clustering and cell definition were performed to identify differentially expressed genes and enrichment pathways. Intercellular communication patterns among major cell types and key signaling pathways were explored by cellular communication analysis.Results:Differential gene expression analysis of FSHD bilateral muscle samples identified 46 functionally differentially expressed genes associated with the disease in different cell types, related to apoptosis, oxidative stress, immune inflammation, and muscle function. Intercellular communication was generally increased in the severe group. Fibro-adipogenic progenitors (FAPs) and macrophages are important signaling sources in the abnormal muscle microenvironment of FSHD and are closely associated with disease progression. There are six unique signaling pathways in the mild group, including bone morphogenetic proteins (BMP), transforming growth factor-β (TGF-β), CXC motif chemokine ligand (CXCL), adhesion G protein-coupled receptor E5 (ADGRE5), interleukin-16 (IL-16), and wingless-type MMTV integration site family (WNT) signaling pathways. These signaling pathways are mainly involved in the interaction between macrophages, FAPs, and adipocytes and may be involved in the regulation of fat deposition and fibrosis changes in the diseased muscle.Conclusions:Single-cell nuclear transcriptome sequencing provides a relatively comprehensive pattern of intercellular communication between key cell types in FSHD, providing an appropriate reference for understanding the intercellular regulatory mechanisms of the FSHD muscle microenvironment.

Loss of ATM activity induces GADD45α-dependent apoptosis in cerebellar granular neurons / 中山大学学报(医学科学版)

objectiveTo explore the specific molecular mechanism of neuronal apoptosis induced by ATM inactivation. MethodsCGNs matured 7 days in vitro were cultured 8 h with 25 K， 5 K or 25 K medium containing ATM-specific inhibitors （Ku55933， 10 µmol/L； Ku60019， 15 µmol/L） for Hoechst stain and apoptosis analysis， or cultured for different lengths of time （2， 4， 8 h） to detect the protein expression levels of ATM， caspase-3 and cleaved caspase-3 by Western blotting. ATM and GADD45α specific siRNA was transfected into C6 cells and CGNs， and its interference efficiency was verified by q-PCR and Western blotting. CGNs matured for 5 days in vitro were transfected with ATM specific siRNA and pCMV-EGFP by calcium phosphate for 48 h， Hoechst staining and apoptosis analysis were performed. CGNs matured for 7 days in vitro were treated with 25 K medium containing ATM specific inhibitors for 8 h， transcriptome sequencing， differential expression gene identification and pathway enrichment analysis were performed. CGNs matured for 5 days in vitro were co-transfected with GADD45α specific siRNA and pCMV-EGFP by calcium phosphate for 48 h， then treated with 5 K or 25 K medium containing 15 µmol/L Ku6 for 8 h. Hoechst staining and apoptosis analysis were performed. ResultsCompared with the 25 K， CGNs nuclear pyknosis rate， cleaved Caspase-3 and ATM protein expression level were increased in the 5 K and ATM-specific inhibitor groups. The mRNA and protein expression levels of ATM and GADD45α were effectively reduced after transfection of ATM and GADD45α specific siRNA in C6 cells and CGNs. Compared with control， CGNs transfected with ATM specific siRNA showed a higher nuclear pyknosis rate. Totally 835 genes were identified to be up-regulated and 848 genes to be down-regulated in the Ku55933 treatment group； 454 genes were identified to be up-regulated and 314 genes to be down-regulated in the Ku6 treatment group； 274 genes were co-up regulated in the Ku5 and Ku60019 treatment groups， while 179 genes were co-down-regulated in the Ku5 and Ku6 treatment groups and the expression of ATM downstream target GADD45α was upregulated. The enrichment results showed that TNF signaling pathway， NF-κB signaling pathway and Apoptosis signaling pathway were significantly enriched. Compared with control， mRNA and protein expression levels of GADD45α were increased in inhibitor treatment and 5 K， while knocking down GADD45α resulted in a decrease in nuclear pyknosis rate in the Ku60019 and 5 K treatment group. ConclusionLoss of ATM activity induces GADD45α-dependent cerebellar granular neuronal apoptosis.

Transcriptomic Effects of LINC01614 on Lung Cancer A549 Cells and Relevance of Drug Resistance / 肿瘤防治研究

Objective To investigate the biological role of LINC01614 in non-small cell lung cancer A549 cells and its drug resistance-related mechanism. Methods The CRISPR/Cas9 technology was used to construct the A549 cell model with knockdown of LINC01614. Transcriptome sequencing was performed on A549 cells knocked down with LINC01614. We validated the transcriptomic differential genes MCAM and ABCC3 at the gene level and MCAM at the protein level, detected the IC50 changes of A549 cells after knockdown of LINC01614 under the effect of different concentrations of cisplatin, and detected the effect of knockdown of LINC01614 on the migration ability of A549 cells. Results Of the 2 713 DEGs after knockdown of LINC01614, a total of 1 626 genes were up-regulated and 1, 087 genes were down-regulated. GO analysis showed that DEGs were associated with intracellular signaling, cell adhesion, and so on. Meanwhile, the KEGG analysis showed that DEGs were associated with the Wnt signaling pathway, TGF-β signaling pathway, and Rap1 signaling pathway. Selection of drug resistance-associated gene ABCC3 from DEGs for validation with MCAM: qRT-PCR results showed that knockdown of LINC01614 significantly down-regulated the expression of MCAM (P<0.05) and upregulated the expression of ABCC3 on A549 cells (P<0.05). After knockdown of LINC01614, the protein expression of MCAM, was significantly decreased in A549 cells (P<0.05); the IC50 of A549 cells to cisplatin was significantly increased (P<0.05); and the scratch healing rate of A549 cells was also significantly decreased (P<0.05). Conclusion LINC01614 may be associated with the proliferation, invasion, and apoptotic pathways of A549 cells. In addition, LINC01614 may exert its migration ability through MCAM and chemoresistance to cisplatin through ABCC3.

Analysis of Tumor Immune Microenvironment Based on Bulk RNA-sequencing Data and Its Application / 肿瘤防治研究

Cancer immunotherapy is one of the most promising biological therapies. The dynamic changes of the immune microenvironment of heterogeneous tumors are critical factors in determining the interaction and therapeutic efficacy between tumor and immune microenvironment. Therefore, quantitative analysis of its constituent cells, related genes, and phenotypes in real time is of great importance. The methods used to analyze tumor immune microenvironment include immunohistochemistry, flow cytometry and so on. With the development of next-generation sequencing technology, transcriptome RNA sequencing data analysis has become one of the important methods to determine the composition of tumor immune microenvironment. This article focuses on the common methods of bulk RNA sequencing data analysis and related research progress.

Mechanism of Huangjing Qianshi Decoction in treatment of prediabetic mice based on transcriptome sequencing / 中国中药杂志

Based on transcriptome sequencing technology, the mouse model of prediabetes treated with Huangjing Qianshi Decoction was sequenced to explore the possible mechanism of treating prediabetes. First of all, transcriptome sequencing was performed on the normal BKS-DB mouse group, the prediabetic model group, and the Huangjing Qianshi Decoction treatment group(treatment group) to obtain differentially expressed genes in the skeletal muscle samples of mice. The serum biochemical indexes were detected in each group to screen out the core genes of Huangjing Qianshi Decoction in prediabetes. Gene Ontology(GO) database and Kyoto Encyclopedia of Genes and Genomes(KEGG) database were used to conduct signaling pathway enrichment analysis of differentially expressed genes, and real-time quantitative polymerase chain reaction(RT-qPCR) was used to verify them. The results showed that the levels of fasting blood glucose(FBG), fasting insulin(FINS), insulin resistance index(HOMA-IR), total cholesterol(TC), triglycerides(TG), and low-density lipoprotein cholesterol(LDL-C) in the mouse model were significantly decreased after treatment with Huangjing Qianshi Decoction. In the results of differential gene screening, there were 1 666 differentially expressed genes in the model group as compared with the normal group, and there were 971 differentially expressed genes in the treatment group as compared with the model group. Among them, interleukin-6(IL-6) and NR3C2 genes, which were closely related to the regulation of insulin resis-tance function, were significantly up-regulated between the model group and the normal group, and vascular endothelial growth factor A(VEGFA) genes were significantly down-regulated between the model group and the normal group. However, the expression results of IL-6, NR3C2, and VEGFA genes were adverse between the treatment group and the model group. GO functional enrichment analysis found that the biological process annotation mainly focused on cell synthesis, cycle, and metabolism; cell component annotation mainly focused on organelles and internal components; and molecular function annotation mainly focused on binding molecular functions. KEGG pathway enrichment analysis found that it involved the protein tyrosine kinase 6(PTK6) pathway, CD28-dependent phosphoinositide 3-kinase/protein kinase B(PI3K/AKT) pathway, p53 pathway, etc. Therefore, Huangjing Qianshi Decoction can improve the state of prediabetes, and the mechanism may be related to cell cycle and apoptosis, PI3K/AKT pathway, p53 pathway, and other biological pathways regulated by IL-6, NR3C2, and VEGFA.

Effect of Flavanomarein from Coreopsis tinctoria Ethanol Extract in Small Intestinal Organoids in Insulin-resistant Mice Based on Transcriptome Sequencing / 中国实验方剂学杂志

ObjectiveTo investigate the effects of flavanomarein on the transcriptome of small intestinal organoids in insulin-resistant mice. MethodFirstly， small intestinal organoids of C57BL/6J and db/db mice were established. Ki-67 and E-cadherin expression was determined by immunofluorescence. Small intestinal organoids were divided into the following three groups： C57BL/6J mouse small intestinal organoids as the normal control group， db/db mouse small intestinal organoids as the model group （IR group）， and db/db mouse small intestinal organoids treated with flavanomarein as the administration group （FM group）. Western blot was used to detect the expression of glucagon-like peptide-1（GLP-1） protein on the small intestinal organoids of the three groups. Finally， transcriptome sequencing was performed on samples from the three groups. ResultOn the 6th day of small intestine organoids culture， a cyclic structure was formed around the lumen， and a small intestine organoids culture model was preliminarily established. Immunofluorescence detection showed that ki-67 and E-cadherin were expressed in small intestinal organoids. Western blot results showed that the expression of GLP-1 protein was increased by flavanomarein. In the results of differential expressed gene （DEG） screening， there were 1 862 DEGs in the IR group as compared with the normal control group， and 2 282 DEGs in the FM group as compared with the IR group. Through protein-protein interaction（PPI） network analysis of the DEGs of the two groups， 10 Hub genes， including Nr1i3， Cyp2c44， Ugt2b1， Gsta1， Gstm2， Ptgs1， Gstm4， Cyp2c38， Cyp4a32， and Gpx3， were obtained. These genes were highly expressed in the normal control group， and their expression was reduced in the IR group. After the intervention of flavanomarein， the expression of the above genes was reversed. ConclusionFlavanomarein may play its role in improving insulin resistance by reversing the expression levels of 10 Hub genes， including Nr1i3， Cyp2c44， Ugt2b1， Gsta1， Gstm2， Ptgs1， Gstm4， Cyp2c38， Cyp4a32， and Gpx3.

Expression pattern of MMP-12 in long-term optic nerve regeneration induced by lens injury / 国际眼科杂志(Guoji Yanke Zazhi)

AIM: To investigate the expression changes of MMP-12 during the long-term axon regeneration induced by the lens injury after the optic nerve clamp trauma in sprague-dawley(SD)rats.METHODS: The optic nerve injury model and lens injury model of SD rats were established, and the 24 experimental animals were divided into control group; lens injury group; optic nerve injury group; lens injury combined with optic nerve injury group, with 6 rats in each group. Reference transcriptome sequencing was used to analyze the expression changes of differentially expressed genes in the injured optic nerve region, and relevant differentially expressed genes with high expression were screened. Quantitative real-time polymerase chain reaction(qRT-PCR)and enzyme-linked immunosorbent assay(ELISA)were used to quantify the expression changes of matrix metalloproteinase-12(MMP-12)in the injured optic nerve region.RESULTS: The Principal Component Analysis of transcriptome sequencing indicated that lens injury combined with optic nerve injury was the principal component of gene expression change. Analysis of gene expression differences showed that the expression of MMP-12 gene was up-regulated in the lens injury combined with optic nerve injury group. The mRNA expression level of MMP-12 in the lens injury combined optic nerve injury group was up-regulated compared with the control group, the optic nerve injury group and the lens injury group at 14d and 21d after successful modeling(P<0.05). At 7, 28d, there was no difference in expression among all groups. The protein expression level of MMP-12 in the lens injury combined with optic nerve injury group was up-regulated compared with the control group and optic nerve injury group at 7, 14 and 21d after successful modeling(P<0.05), and it was up-regulated in the lens injury group combined with optic nerve injury group compared with optic nerve injury group at 21d(P<0.05). At 28d, there was no difference in expression among all groups.CONCLUSION: The up-regulated expression of MMP-12 may be involved in the long-term regeneration of the optic nerve after lens injury.

Mechanism of Huashi Baidu Prescription in Treating Influenza Viral Pneumonia Based on Transcriptome / 中国实验方剂学杂志

ObjectiveTo clarify the therapeutic effect of Huashi Baidu prescription on pneumonia in mice caused by influenza A （H1N1） virus and explore its mechanism based on the transcriptome. MethodA mouse influenza viral pneumonia model was built by intranasal infection with influenza A virus， and mice were continuously administered the drug for five days， so as to investigate the general condition， lung index， viral load， pathological morphology of lung tissue， survival time， and prolongation rate of survival time of mice and clarify the therapeutic effect of Huashi Baidu prescription on influenza viral pneumonia. Transcriptome technology was used to detect the differentially expressed genes in the lung tissue of mice in the model group and the normal group， as well as the Huashi Baidu prescription group and the model group， and the potential core target of the Huashi Baidu prescription for the treatment of influenza viral pneumonia was screened. Real-time fluorescence quantitative polymerase chain reaction （Real-time PCR） was used to verify the effect of Huashi Baidu prescription on the mRNA expression level of core target genes. ResultCompared with the normal group， the lung index and viral load in the lung tissue of the model group were significantly increased （P<0.05， P<0.01）. Compared with the model group， the high-dose group of Huashi Baidu prescription significantly prolonged the survival time of mice infected with influenza A virus （P<0.05） and significantly reduced the lung index value of mice （P<0.05） and the viral load of lung tissue. The high-dose， medium-dose， and low-dose groups of Huashi Baidu prescription could significantly reduce lung tissue inflammation， blood stasis， swelling， and other pathological changes in mice （P<0.05， P<0.01）. Transcriptome analysis of lung tissue showed that core genes were mainly enriched in the nuclear transcription factor-κB （NF-κB） signaling pathway， interleukin-17 （IL-17） signaling pathway， cytokine-cytokine receptor interaction， and other pathways after the intervention of Huashi Baidu prescription. TRAF6， NFKBIA， CCL2， CCL7， and CXCL2 were the top five node genes with combined score values. Real-time PCR validation showed that Huashi Baidu prescription significantly downregulated the mRNA expression of key genes TRAF6 and NFKBIA in the NF-κB signaling pathway， as well as chemokines CCL2， CCL7， and CXCL2 （P<0.05， P<0.01）. ConclusionHuashi Baidu prescription has a therapeutic effect on influenza viral pneumonia， possibly by inhibiting the expression of key nodes TRAF6 and NFKBIA in the NF-κB signaling pathway and that of chemokines CCL2， CCL7， and CXCL2， reducing the recruitment of inflammatory cells and viral load， and exerting anti-influenza viral pneumonia effects.

Construction of CeRNA network of differential LncRNAs in THP-1 cells infected by Dengue virus type Ⅲ and its antibody dependent enhancement / 中国生物制品学杂志

@#ObjectiveTo establish models of Dengue virus type Ⅲ（DENV-3,DV-3）infection and antibody dependent enhancement(ADE)infection at the acute monocytic leukemia cells(THP-1),investigate the differential expression of long non-coding RNAs(LncRNAs),map the competitive endogenous RNA(CeRNA)regulatory network and predict the translation function of LncRNAs.MethodsThe culture supernatant was harvested 6 d after C6/36 cells were infected with DENV-3,the virus titer was determined by CCID50,and the type and full-length genome amplification were identified by PCR;The DENV-3 standard plasmid was amplified,identified by PCR,and the standard curve was drawn;THP-1 cells were divided into negative control group(THP-1),direct infection group(DV-3),ADE group and blank control group[1640(-)]. After 48 h of infection,the total RNA was extracted and the copy number of intracellular virus nucleic acid was measured;Through the whole transcriptome sequencing technology,the CeRNA regulatory network was constructed for the top five up-regulated and down-regulated LncRNAs in THP-1 vs DENV3,THP-1 vs ADE,DENV3 vs ADE groups,and the functions of their coding proteins were analyzed.ResultsC6/36 cells infected with DENV-3 for 3 d showed obvious cell fusion,vacuoles and abscission;The virus had a titer of about 1. 0 × 104. 64PFU/mL and was identified as DENV-3 by PCR specific primers,of which the complete gene sequence was obtained;The number of viral nucleic acid copies in ADE group was significantly higher than those in DV-3 group and blank control group;In THP-1 vs DENV-3,the expression of cytohesin interacting protein(CYTIP)was predicted to be up-regulated;In THP-1 vs ADE,the expression of kinesin family5A(KIF5A)was predicted to be down-regulated;In DENV-3 vs ADE,the expression of cluster differentiation antigen 9(CD9)and insulin like growth factor 2(IGF2)was predicted to be up-regulated. All of these differential LncRNAs had open reading frames(ORFs). Except Lnc-SH3BP1 and Lnc-RPL41,all of the other LncRNAs had internal ribosome binding site(IRES).ConclusionIn DENV-3 infection of THP-1 cells and ADE infection mediated by DENV-3,the expression of LncRNAs has changed significantly,and may regulate the process of infection through a variety of biological functions,which is helpful for a deeper understanding of the mechanism of ADE infection.

Transcriptome sequencing of transgelin-2 inhibiting high glucose induced microglia inflammation / 中华眼底病杂志

Objective:To analyze the change of differential genes and signaling pathways in high glucose induced BV2 cells, and to explore the mechanism of transgelin-2 (TAGLN2) regulating cellular inflammatory response and metabolic process.Methods:An experimental study. The cultured BV2 cells were divided into mannitol treatment (Man) group, glucose treatment (Glu) group, overexpression control Glu treatment (Con) group, overexpression TAGLN2 Glu treatment group, silence control Glu treatment (shCon Glu) group, and silence TAGLN2 Glu treatment (shTAGLN2 Glu) group. Cells in the Man group were cultured in modified Eagle high glucose medium (DMEM) containing 25 mmol/L mannitol and 25 mmol/L glucose, cells in other groups (Glu group, Con Glu group, TAGLN2 Glu group, shCon Glu group and shTAGLN2 Glu group) were cultured in DMEM medium containing 50 mmol/L glucose. After 24 hours of cells culture, transcriptome sequencing of cells in each group were performed using high-throughput sequencing technology, and significantly differentially expressed genes (DEG) were screened. |log 2 (fold change)|≥1 and P≤0.05 were adopted as criteria to screen for DEG. Gene Ontology (GO) function enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis and protein-protein interaction network analysis were performed. Real-time polymerase chain reaction (RT-PCR) was used to detect the relative expression level of DEG mRNA. The data between groups were compared by independent sample t-test. Results:When compared with Man group, a total of 517 differentially expressed genes were screened in Glu group, which including 277 up-regulated genes and 240 down-regulated genes. KEGG pathway enrichment analysis showed that the up-regulated genes were significantly enriched in immune system processes such as nuclear factor (NF)-κB signal pathway, Jak-signal transducers and activators of transcription (STAT) signal pathway, while down-regulated genes were significantly enriched in glycosaminoglycan degradation and glyceride metabolic pathway. Compared with Con Glu group, a total of 480 DEG were screened in TAGLN2 Glu group, among which 147 up-regulated and 333 down-regulated genes were detected. Up-regulated genes were significantly enriched in the metabolic processes of fatty acid, glyceride and pyruvate, while down-regulated genes were significantly enriched in immune system processes such as NF-κB signal pathway, Jak-STAT signal pathway and tumor necrosis factor (TNF) signal pathway. Compared with shCon Glu group, a total of 582 DEG were screened in shTAGLN2 Glu group, among which 423 up-regulated and 159 down-regulated genes were detected. Up-regulated DEG were significantly enriched in immune system processes such as TNF signal pathway and chemokine signal pathway, while down-regulated DEG were significantly enriched in pattern recognition receptor signal pathway. RT-PCR results showed that the relative expression levels of DEG mRNA Card11 ( t=13.530), Icos ( t=3.482), Chst3 ( t=6.949), Kynu ( t=5.399), interleukin (IL)-1β ( t=2.960), TNF-α ( t=5.800), IL-6 ( t=3.130), interferon-γ ( t=7.690) and IL-17 ( t=6.530) in the TAGLN2 Glu treatment group were decreased significantly compared with Con Glu group, and the difference was statistically significant. Conclusion:TAGLN2 can inhibit glucose induced microglia inflammation by NF-κB and Jak-STAT signaling pathways, Card11, Icos, Chst3 and Kynu play an important role in the anti-inflammatory process of TAGLN2.

Specific expression of FABP5 in type II alveolar epithelial cells in silicosis pulmonary fibrosis based on single-cell transcriptome sequencing / 环境与职业医学

Background Silicosis is caused by long-term inhalation of large amounts of free silica (SiO2) particles, and exploring its mechanism can provide new directions for the treatment of silicosis fibrosis. Objective To investigate the expression and role of fatty acid binding protein 5 (FABP5) in a silica-induced silicosis model. Methods In combination with the results of single-cell transcriptome sequencing, the expression pattern of FABP5 in mouse alveolar epithelial cells was explored by bioinformatic analysis, and the distribution pattern of fabp5 was detected by spatial transcriptomics. An in vivo model of silicosis was established by intratracheal injection with SiO2 into mice and four groups were set up: normal saline (NS) 7 d group, NS 56 d group, SiO2 7 d group, and SiO2 56 d group. An in vitro model of silicosis was established in SiO2-treated mouse lung epithelial cell line (MLE-12). At the whole animal level, the marker of epithelial cells (E-Cad) and the protein level of FABP5 were detected by tissue immunofluorescence assay; in vitro, the changes of fabp5 mRNA expression and protein level in MLE-12. Results The results of single-cell transcriptome sequencing and spatial transcriptome sequencing showed that the mRNA expression of fabp5 was upregulated in type II alveolar epithelial cells in the focal area of silicosis in mice, accompanied by elevated tissue immunofluorescent protein levels, and there was co-localization of E-CAD. Meanwhile, SiO2 stimulation induced a 1.58-fold increase in fabp5 mRNA expression and a 2-fold increase in protein levels in MLE-12 cells, with significant differences (P<0.05). Conclusion The protein level of FABP5 is increased in alveolar epithelial cells in a pulmonary fibrosis model, suggesting that FABP5 may be involved in the pathological process of epithelial cells in pulmonary fibrosis.

Functions of key genes involved in TGF-β/Smad signaling pathway in progression of pulmonary fibrosis / 环境与职业医学

Background Although transforming growth factor-β (TGF-β)/Smad signaling pathway is important in regulating the occurrence and development of pulmonary fibrosis, the pathogenesis of pulmonary fibrosis remains elusive. Objective To explore the functions of genes associated with TGF-β/Smad signaling pathway in the progression of pulmonary fibrosis. Methods A NIH-3T3 fibroblast model induced by TGF-β1 was established. The experiment samples were divided into a control group and a TGF-β1 treatment group. The control group was exposed to normal saline, while the TGF-β1 treatment group was exposed to 10 ng·mL−1 TGF-β1 for 12 h. The RNAs of the two groups were extracted, sequenced, and analyzed by bioinformatics methods to identify seven key genes in TGF-β pathway, including Dcn, Smad3, Smad7, Fbn1, Thbs1, TGF-β1, and TGF-β3. The gene expression levels of five markers [Collagen1α1, Collagen1α2, α-smooth muscle actin (α-SMA), TGF-β1, and TGF-β3] and the seven key genes were detected by quantitative real-time PCR (qRT-PCR). The proteins of the two groups were extracted. The important marker protein expression levels of Smad3, the phosphorylation of Smad3 (P-Smad3), and α-SMA were detected by Western blotting. At the same time, 30 healthy SPF-grade C57BL/6 mice were randomly divided into three groups, with 10 mice in each group: a control group, a SiO2 inhalation exposure group for 28 d (10 mice), and a SiO2 inhalation exposure group for 56 d (10 mice). The mice in the two treatment groups were exposed to a natural SiO2 environment for 4 h per day with a 10-min pause for breathing fresh air at 2 h intervals. The lung tissues of the mice were taken after execution. The changes of pulmonary fibrosis were detected by Masson staining, and mRNAs and proteins were extracted to detect the expression of the above key genes and proteins. Results The expression levels of the five marker genes Collagen1α1, Collagen1α2, α-SMA, TGF-β1, and TGF-β3 were significantly increased in the TGF-β1-induced NIH-3T3 fibroblasts than those in the control group (P < 0.01); the expression levels of P-Smad3 and α-SMA proteins increased significantly (P < 0.01); the expression results of the seven key genes screened in the TGF pathway were that Dcn and Smad3 were obviously down-regulated (P < 0.01), and Smad7, Fbn1, Thbs1, TGF-β1, and TGF-β3 were obviously up-regulated (P < 0.01). The changes in gene expression levels of the transcriptome sequencing showed the same trend. The results of Masson staining showed that the content of collagen fibers in the lung tissues also increased in the SiO2 inhalation exposure groups over time. In the mouse experiment, five marker genes were obviously up-regulated compared with the control group (P < 0.01); no obvious change was found in the expression of Smad3 protein, and the expression levels of P-Smad3 and α-SMA were obviously higher in the SiO2 exposure groups than those in the control group (P < 0.01); the expression levels of Dcn and Smad3 showed a down-regulated trend, while the expression levels of Smad7, Fbn1, Thbs1, TGF-β1, and TGF-β3 showed an up-regulated trend with the increase of SiO2 inhalation exposure days (P < 0.01). The expression levels of the above five marker genes, three important marker proteins, and seven key genes were consistent with the expression trends of TGF-β1-induced NIH-3T3 fibroblasts. Conclusion The expression levels of pulmonary fibrosis-related marker genes and proteins change significantly in TGF-β1-induced fibroblast cells, and the lung tissues of mice under natural SiO2 inhalation exposure has obvious fibrosis characteristics. Seven genes (Dcn, Smad3, Smad7, Fbn1, Thbs1, TGF-β1, and TGF-β3) may be involved in the regulation of pulmonary fibrosis by the TGF-β/Smad signaling pathway.

Analysis ofthe Mechanisms of ZNF185in Esophageal Squamous Cell Carcinoma Based on Transcriptome Sequencing / 华中科技大学学报(医学版)

Objective To explorethe mechanism and down stream signaling pathways of ZNF185 in esophageal squamous cell carcinoma by transcriptome sequencing and bioinformatics analysis. Methods The expression of ZNF185 was knocked down by conducting shRNA transfection of esophageal squamous carcinoma cell lines,and the changes of mRNA expression were detected by transcriptome sequencing technology. The differential genes were screened based on sequencing data,and the differential gene-related functions and signaling pathways were analyzed by GO functional analysis and KEGG enrichment analysis. The interactions of differential proteins were analyzed by STRING,and the related protein-protein interaction(PPI) network was constructed to analyze the keynode proteins. Results There were 3859 differential genes between shZNF185 group and the control group,including 2112 genes with elevated expressionand 1747 genes with decreased expression. RNA splicing cell-matrix junction,local adhesion were significantly enriched in GO analysis. KEGG enrichment analysis indicated that ZNF185 might be related to HIF-1 signaling pathwayand p53 signaling pathway,and key proteins such as BMP2 and OAS2 were identified in PPI network. Conclusion The possible signaling pathways involved in the regulation of esophageal squamous carcinoma by ZNF185 werescreened out by transcriptome sequencing and bioinformatic analysis, which provide a newidea to elucidate the pathogenesis associated with esophageal squamous carcinoma.

Whole-transcriptome sequencing analysis of placental differential miRNA expression profile in Down syndrome / 南方医科大学学报

OBJECTIVE@#To identify new biomarkers and molecular pathogenesis of Down syndrome (DS) by analyzing differentially expressed miRNAs in the placentas and their biological pathways.@*METHODS@#Whole transcriptome sequencing was used to identify the differentially expressed miRNAs in DS (n=3) and normal placental samples (n=3) diagnosed by prenatal diagnosis. The target genes were predicted using miRWalk, Targetscan and miRDB, and GO and KEGG pathway analyses were performed for gene enrichment studies.@*RESULTS@#We identified a total of 82 differentially expressed miRNAs in the placental tissues of DS, including 29 up-regulated miRNAs (fold change ≥2, P < 0.05) and 15 down-regulated miRNAs (fold change ≥2, P < 0.05), among which 10 miRNAs with relatively high expression abundance were selected for further analysis, including 4 up-regulated and 6 down-regulated miRNAs. These selected miRNAs shared the common target genes BTBD3 and AUTS2, both of which were associated with neurodevelopment. GO analysis showed that the target genes of the selected miRNAs were mainly enriched in protein binding, hydrolytic enzymes, metal ion binding protein combining, transferase activity, nucleotide, cytoplasmic constituents, nucleus composition, transcriptional regulation, RNA metabolism regulation, DNA-dependent RNA polymerase Ⅱ promoter transcriptional regulation, eye development, and sensory organ development. KEGG enrichment analysis showed that the target genes of these differentially expressed miRNAs were involved in the signaling pathways including tumor-related signaling pathway, PI3K-Akt signaling pathway, Ras signaling pathway, Rap1 signaling pathway, cytoskeletal regulatory signaling pathway, purine metabolization-related signaling pathway and P53 signaling pathway.@*CONCLUSION@#The differentially expressed miRNAs may play important roles in placental damage and pregnancy pathology in DS and provide clues for the prevention and treatment of mental retardation-related diseases.

Application of single-cell RNA-sequencing in age-related macular degeneration / 国际眼科杂志(Guoji Yanke Zazhi)

@#Age-related macular degeneration(ARMD)is the most common eye disease that can cause irreversible loss of central vision in the elderly population. Since the complexity of the pathogenesis in ARMD, the underlying mechanism remains uncovered,and treatment limited. Conventional bulk transcriptome sequencing strategies can only provide the average gene profile in the dominant cells, while single-cell RNA-sequencing(scRNA-seq)is able to reveal the mRNA transcriptome at a single-cell level. The ScRNA-seq has been applied to discover novel cell subtypes, reveal intercellular heterogeneity, and unveil the process of cell differentiation.In this paper, we reviewed the technical principle of scRNA-seq and its application in retinal,choroid development and ARMD research, raised the defection of scRNA-seq and trended in emerging technologies,provided the new insights and perspectives for the in-depth study of retinal and choroid physiology, pathophysiology, and pathogenesis of ARMD diseases. It is hoped to provide theoretical foundation for the targeted therapy of ARMD.

Whole transcriptome analysis and critical gene regulatory network analysis during Schistosoma japonicum infection and praziquantel treatment in mice / 中国血吸虫病防治杂志

Objective To investigate long non-coding RNA (lncRNA)-microRNA (miRNA)-messenger RNA (mRNA) interactions and identify the critical gene regulatory network during Schistosoma japonicum infections and praziquantel treatment using whole transcriptome sequencing. Methods A total of 110 male C57BL/6 mice were randomly divided into the control group, the infection group and the treatment group. Mice in the infection treatment and the control group were infected with S. japonicum cercariae via the abdomen, and liver specimens were sampled from 10 mice 3, 6, 8 weeks post-infection. Praziquantel treatment was given to mice in the treatment group 8 weeks post-infection, and liver specimens were sampled from 10 mice 2, 4, 6, 8, 10 weeks post-treatment. Total RNA was isolated from mouse liver specimens, and the transcriptome library was constructed for highthroughput whole transcriptome sequencing. The significant differentially expressed genes were subjected to functional annotations, Gene Ontology (GO) terms enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Correlation analysis of liver specimens was performed using R Corrplot and Himsc functions, and the lncRNAmiRNA-mRNA interaction network analysis was performed using R MixOmics and Himsc functions. Results There were 1 176 differentially expressed miRNAs, 5 270 differentially expressed mRNAs, and 2 682 differentially expressed lncRNAs between the infection group and the control group, 1 289 differentially expressed miRNAs, 7 differentially expressed mRNAs, and 69 differentially expressed lncRNAs between the treatment group and the infection group, and 1 210 differentially expressed miRNAs, 4 456 differentially expressed mRNAs, and 2 016 differentially expressed lncRNAs between the treatment group and the control group. Correlation analysis showed a higher correlation of gene expression between the treatment group and the control group. Principal component analysis showed obvious separate clustering between the infection group and the treatment group. The differentially expressed genes with significant relevance were significantly enriched in 24 GO terms, including arachidonic acid metabolic process, xenobiotic catabolic process, unsaturated fatty acid metabolic process, xenobiotic metabolic process, long-chain fatty acid metabolic process, and 8 KEGG metabolic pathways, including cholesterol metabolism, tyrosine metabolism, linoleic acid metabolism, retinol metabolism, and steroid hormone biometabolism. Conclusions There were 23 mRNAs including Cyp2b9 and 14 lncRNAs including Rmrpr in the core position of the gene regulatory network, which may play a critical role in S. japonicum infections and praziquantel treatment, and 9 miRNAs including miR-8105 may serve as potential molecular markers for diagnosis of S. japonicum infections.

Acute leukemia with NUP98-DDX10 fusion gene-positive: report of 2 cases and review of literature / 白血病·淋巴瘤

Objective:To investigate the clinical characteristics, diagnosis, treatment and prognosis of acute leukemia (AL) with NUP98-DDX10 fusion gene-positive.Methods:The clinical data of 2 AL patients with NUP98-DDX10 fusion gene-positive who admitted to Blood Diseases Hospital, Chinese Academy of Medical Sciences in April 2020 and February 2021, respectively were retrospectively analyzed. Transcriptome gene sequencing was used to detect fusion gene, and the fusion gene fragment was amplified by using reverse transcription polymerase chain reaction (RT-PCR), and Sanger sequencing was used to clarify sequences. The clinical and experimental indicators characteristics were analyzed and the relevant literatures were reviewed.Results:According to the clinical diagnosis, 1 patient was diagnosed as acute myeloid leukemia M 5 (AML-M 5) and 1 patient was diagnosed as acute leukemia of ambiguous lineage, not otherwise specified (ALAL-NOS). The AML-M 5 patient presented with severe coagulation abnormalities, and fulfilled the diagnostic criteria for diffuse intravascular coagulation (DIC) at the initial visit. Transcriptome sequencing of 2 patients showed NUP98-DDX10 fusion gene- positive. RT-PCR confirmed that sequencing results identified 2 different splice fusion modes: one was NUP98 exon 14 fused with DDX10 exon 7(usually called "type Ⅱ"), the other was NUP98 exon 14 fused with DDX10 exon 13, which was never reported and named as "type Ⅳ". From 1997 to 2018, a total of 16 cases with NUP98-DDX10 related hematologic neoplasms were reported in the literature. A summary analysis of 16 cases added with 2 patients in our center included 13 males and 5 females with median age 31.5 years (0.08-61 years). The median overall survival was 12 months (1-46 months). Conclusions:A novel fusion gene NUP98-DDX10 transcriptome is identified in ALAL-NOS patient. Hematological malignancies with NUP98-DDX10 are very rare. They respond poorly to conventional treatment and require allogeneic hematopoietic stem cell transplantation (allo-HSCT) to improve the prognosis.